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Abstract 

This paper studies the Yang-Lee edge singularity of 2-dimensional 
(2D) Ising model based on a quantum spin chain and transfer matrix 
measurements on the cylinder. Based on finite-size scaling, the low-lying 
excitation spectrum is found at the Yang-Lee edge singularity. Based 
on transfer matrix techniques, the single structure constant is evaluated 
at the Yang-Lee edge singularity. The results of both types of measure- 
ments are found to be fully consistent with the predictions for the (A4, A\) 
minimal conformal field theory, which was previously identified with this 
critical point. 

1 1. Introduction 

In 1978, Fisher [1] proposed that Yang-Lee edge singularities [2, 3] are critical 
points. Later, Cardy [4] argued that the Yang-Lee edge singularity of the 2D 
Ising model should be identified with the (A^^Ai) minimal conformal field the- 
ory (CFT) [5, 6] of the ADE classification [7]. Cardy's identification provides 
CFT predictions for this Yang-Lee edge singularity. 

This article tests different predictions coming from Cardy's identification. 

In section 2, we provide measurements of the low- lying excitation spectrum 
at Yang-Lee edge singularity of the 2D Ising model. The measured low-lying 
excitation spectrum is also compared with predictions from Cardy's identifica- 
tion of the (A4, Ai) minimal CFT with this Yang-Lee edge singularity of the 2D 
Ising model [4, 8, 9]. 
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Cardy's identification also determines the forms of 2-point and 3-point cor- 
relations. In particular, these correlations define universal amplitudes, which 
are known as structure constants [5, 10]. Such predictions are an important ad- 
vance that CFT brought to the understanding of critical points of 2D statistical 
models. No tests of such predictions have been performed for critical points 
associated with non-unitary CFTs. 

In section 3, we provide a measurement of the universal amplitude associ- 
ated with the Yang-Lee edge singularity of the 2D Ising model. The measured 
amplitude is also compared with the prediction from Cardy's identification of 
the (Ai,j4i) minimal CFT with this Yang-Lee edge singularity. 

2 2. Excitation Spectrum at the Yang-Lee Edge 
Singularity of the 2D Ising model 

The 2D Ising model in an imaginary external magnetic field is associated with 
a quantum spin chain whose Hamiltonian, Hi S i ng , on an iV-site chain, is given 
by [15]: 

N 

Hi s i n g = - ^2{t(J z (n)<T z (n + 1) + iBa z (n) + <r x (n)}. (1) 

n=l 

In Eq. (1), <j x (n) and cr z (n) are Pauli spin matrices at the site n, parameter 
n t n is a positive coupling for a ferromagnetic spin-spin interaction, and iB is a 
purely imaginary external magnetic field. In Eq. (1), the last term produces 
inter-row single spin flips in the associated 2D transfer matrix [16, 17]. 

Below, the phenomenological renormalization group (PRG) is used to deter- 
mine critical values of imaginary magnetic field, iB Y L,{N), for various lengths, 
N, of the chain. For imaginary magnetic fields, the PRG equation requires 
that: [11, 8] 

[N - l}m{B Y L(N),N - 1) = [N]m(B YL (N),N). (2) 

In Eq. (2), m(B,N) = [Ei(B,N) -E (B,N)] where E (B,N) and Ei(B,N) 
are energies for the ground state "0" and the first excited state "1" on a chain 
of length N. Below, m(B, N) is referred to as Gap(B, N) or more simply as 
Gap(N). At these B Y l(N)'s, the Ising quantum spin chain exhibits the finite- 
size scaling behavior of the Yang-Lee edge singularity. In particular, if the 
ByLiNys converge to a nonzero value as N — > oo, that value will be the critical 
point for the Yang-Lee edge singularity of the 2D Ising spin model. 

At these Byl{N)'s, excitation energies should scale. In particular, CFT 
predicts how these energies will scale with the length, N, of the chain. For an 
excited energy eigenstate " i" of the quantum spin chain, an excitation energy, 
Ei{N) - Eq(N), will scale as: [14] 

8W - W .» ^M . (3) 
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In Eq. (3), Aj and Aj are left and right conformal dimensions of conformal 
field " i" , and A and A are conformal dimensions of the primary field having the 
lowest "negative" scaling dimension in the relevant non- unitary CFT. In Eq. 
(3), the constant £ is non-universal, e.g., depending on the normalization of the 
Hamiltonian . 

In minimal CFTs, the modular invariant forms of the partition functions 
[7, 12] determine the low-lying excitation spectrum and the central charges. For 
the (A4, A\) minimal CFT, Table 1 gives the energies of the low-lying excitations 
and degeneracies thereof as obtained from the associated partition function. 
Table 1 provides normalized excitation energies, which are ratios. For a state "i" , 



CFT 


(A 4 ,A 1 ) 




Normalized Energies 





1 


2.5 


5.0 


6.0 


7.5 


Degeneracy 


1 


1 


2 


3 


2 


4 



Table 1: Lowest excitations of (Aj,^4i) CFT 



the normalized excitation energy is the ratio is the excitation energy of the state 
"i" over the excitation energy of the lowest excited state "1". Here, excitation 
energies are with respect to the ground state. The normalized excitation energies 
of Table 1 do not depend on non-universal constants such as £. 

The critical magnetic fields, Byl(N), were obtained by solving the PRC 
eq. (2) for chains of different lengths. For these solutions, state energies were 
obtained by using the Lanczos algorithm for Hj S i ng of eq. (1). Table 2 shows 
critical fields, i.e., Byl{N)'s, ground state energies, and lowest excitation ener- 
gies, i.e., Gap(TV)'s. These measurements were obtained for Ising quantum spin 
chains in which the coupling, t, is 0.1. Table 2 shows that NxGap(N) scales to 
a constant as TV — > oo as expected from the PRC. 

The critical magnetic field values, i.e., the Byl's, were used to find the low- 
lying excitation spectra of Ising quantum spin chains of various lengths. Table 
3 provides measured spectra including both energies and degeneracies. Here, 
excitation energies are also normalized by dividing by the lowest excitation 
energy, i.e., as already described to remove any dependence on the non-universal 
constant £. 

Figures 1-4 plot the measured excitation energies of the states A - F as a 
function of the inverse of the length of the Ising quantum spin chain. 

A visual inspection of Figures 1-4 shows that the type A, B, C, D, E, and 
F states form four distinct sets A, B & C, D, and E & F. Within each set, the 
states have energies that approach the same value as 1/TV — ► 0. The excitation 
energies of the states of sets A, B & C, D, and E & F approach about 2.45, 5.0, 
6.03, and 7.6, respectively, as 1/TV — > 0. A BST analysis shows that excitation 
energies of the type A, B, C, D, E, and F states scale to 2.4995(5), 5.005(1), 
5.003(3), 5.99(1), 7.54(8), and 7.60(7), respectively, in this limit. These PRC 
measurements of the low-lying excitation energies and degeneracies agree well 
with the predictions for the (A±, A±) CFT as in Table 1. 

1 iJ is the "sound velocity" in the dispersion relation of the critical Hamiltonian 
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Figure 1: Energies of the type A states as a function of 1/N. 
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Figure 2: Energies of type B states (squares) and type C states (circles) as a 
function of 1/iV. 
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Figure 3: Energies of type D states as a function of 1/iV. 
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Figure 4: Energies of type E states (squares) and F states (circles) as a function 
of 1/N. 
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l\Tn m nPT" 

of Sites 


Byl(N) 


Energy of ground state 


Gap(N) 


N x Gap(N) 


3 


.2459180i 


-2.8811043 


.8103423 


2.4310 


4 


.2384127i 


-3.8028211 


.6629112 


2.6516 


5 


.2352339i 


-4.7341982 


.5613016 


2.8065 


6 


.2337637i 


-5.6688215 


.4858628 


2.9152 


7 


.2330279i 


-6.6048003 


.4275400 


2.9928 


8 


.2326347i 


-7.5414746 


.3811698 


3.0494 


9 


.2324118i 


-8.4785910 


.3435105 


3.0916 


10 


.2322793i 


-9.4160213 


.3123765 


3.1237 


11 


.2321972i 


-10.353696 


.286250 


3.1488 


12 


.2321442i 


-11.291568 


.264041 


3.1685 


00 


.23193i 


— oo 


0.0 


3.2840 



Table 2: Measurements of Byl(N), Ground state energy, Gap, and NxGap for 
various chain lengths, N 



State /[Degeneracy] 


6 


7 


8 


9 


10 


11 


12 


A/ [2] 


2.68432 


2.64386 


2.61415 


2.59207 


2.57540 


2.56260 


2.55253 


B/[l] 


4.18193 


4.27896 


4.36713 


4.44474 


4.51197 


4.56977 


4.61912 


C/[2] 


4.51738 


4.63236 


4.70368 


4.75182 


4.78652 


4.81281 


4.83329 


D/[2] 


5.85889 


5.89208 


5.91240 


5.92644 


5.93703 


5.94544 


5.95210 


E/[2] 




5.68559 


6.03104 


6.27270 


6.45018 


6.58573 


6.69223 


F/[2] 




6.24252 


6.35798 


6.46344 


6.55966 


6.64694 


6.72535 



Table 3: Normalized excitation energies and degeneracies of lowest excited states 
A - F for Ising quantum spin chains with 6 to 12 sites. 



3 3. Structure Constant at the Yang-Lee Edge 
Singularity 

The non-unitary (Ai,Ai) minimal CFT has one primary field <p{z, z) with left 
and right conformal weights -1/5 and scaling dimension x of -2/5 [4]. For this 
field, (j)(z, z), 2-point and 3-point correlations have the forms: 

G lf>(t> (z 1 ,z 1 ,z 2 ,z 2 ) = \{zi - z 2 )\ 4/5 , (4) 

and 

G^izi, zi, z 2 , z 2 , z 3 ,z 3 ) = C\(zi - z 2 )(z 2 - z 3 )(z 3 - zi)| 2/5 . (5) 

Cardy showed that the structure constant, C, of the non-unitary (A±, Ai) min- 
imal CFT is given by: [4] 

/ [r(6/5)Ml/5)r(2/5) 
V r(3/5)[T(4/5)]3 1 > 
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Below, numerical measurements at the Yang-Lee edge singularity are presented 
for this CFT prediction. The numerical measurements were made for the 2D 
Ising model, i.e., rather than for a spin chain. The 2D Ising model has a Hamil- 
tonian H, given by: 

M N 

H = — 2^ + + hSij]. (7) 

j=l i=l 

The spin-spin coupling J is positive. In this model, the Yang-Lee edge singu- 
larity occurs above the critical temperature for a purely imaginary values of the 
magnetic field, h, i.e., h = iB with B real [2, 3]. Below, the spin correlations 
were measured at a temperature, T, for which J/ksT = 0.1 

The transfer matrix was used to measure correlation 2-spin and 3-spin cor- 
relations on torii of length, M, and diameters, N. In these evaluations, M was 
much larger than N, i.e., M = 512 and N = 3 - 8, so that correlations had dis- 
tance behaviors for infinitely long cylinders at field separations small compared 
to M. 

Finite-size scaling enabled the extraction of physical properties in the ther- 
modynamic limit [13]. In particular, the spin correlations were measured at 
purely imaginary magnetic field values, h(N) = iByl{N). Each value , Byl(N), 
satisfied the phenomenological renormalization group (PRG) equation for infi- 
nite cylinders of diameters (N — 1) and N: 

S(iB YL (N),N-l) = giB YL (N),N) 

TV- 1 N [ ' 

In the PRG equation, £(iB, N) is the spin-spin correlation length on the infinite 
cylinder of diameter N at the magnetic field iB 2 . 

On a cylinder of width N, CFT predicts that correlations depend exponen- 
tially on distances between fields when said distances are large compared the 
cylinder's diameter, N [14]. When \yi — J/2I >> N, the 2-point correlation of 
fields of scaling dimension, x, has the form exp(— 2-Kx(yi — 1)2) /N) where y\ 
and J/2 are the positions of the fields along the axis of the infinite cylinder. For 
the 3-point correlation, the exponential behavior on the distances between the 
fields of the correlation is also determined by the scaling dimensions of the fields 
therein 

At the Yang-Lee edge singularity, amplitudes of 2-spin and 3-spin correla- 
tions, i.e., A ss and A sss , were used to evaluate the 3-spin structure constant. 
The values of the 3-spin structure constant, C(N), were obtained from the re- 
lation: 

A SSS ^B YL {N)) 
° liVj [A sa {iB YL {N)]W W 

In the above equation, A ss (iB Y ^ (N)) and A sss (iB Y ^ (N)), are amplitudes of the 
respective 2-spin and 3-spin correlations at PRG values for the magnetic field. 



2 £ is measured by the first inverse gap (see [16]) 
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PRG measurements of the correlation length £(N) also provide a measurement 
of the conformal dimension, x, of the spin field, i.e., x(N) = N/[27r£(iBYL {N))} ■ 
Scaling behaviors of these quantities with the cylinder's width, N, were used to 
obtain the values of the quantities as N — > oo. 

Table 4 provides our transfer matrix measurements 3 for M = 512 and 
J/k B T = 0.1. 



N 


Byl(N) 


x(N) 


\C(N)\ 


3 


0.184802 


0.353929 


1.80838 


4 


0.183348 


0.376870 


1.83711 


5 


0.183064 


0.385748 


1.85736 


6 


0.182982 


0.390108 


1.87054 


7 


0.182951 


0.392693 


1.87937 


8 


0.182946 


0.392911 


1.88633 


00 




0.398(2) 


1.923(13) 


CFT 




0.4 


1.9113 



Table 4: PRG Measured values of conformal dimension and structure constant. 

In Table 4, the oo line extrapolates the measured values to the thermody- 
namic limit, i.e., N = oo. The extrapolated values were obtained from fits of 
the measured x(N) y s and \C(N)\'s to functions of form f(N) = /(oo) + fi N~ a , 
i.e., leading finite-size scaling forms. 

In Table 4, the last line gives the predictions for x and \C\ from the (Aj, A) 
non-unitary minimal CFT model. 

Figure 5 shows measurements of the structure constant |C(7V)| and a best 
fit (line), which accounts for a correction in TV -12 . In Figure 5, the black 
squares are the measured C(N)'s, the empty square is the value of C(oo) from 
the best fit, and the black circle is the CFT prediction. Our finite-size scaling 
measurements produce a value for the structure constant that again agrees well 
with the prediction of the (A4, A) minimal CFT. 



3 For N = 3 - 7, B Y l{N)'s were evaluated from the PRG equation. For N = 8, Byl(S)' was 
estimated from the B Y £ y (N) , s for N = 3 - 7 by assuming a leading finite-size scaling behavior. 
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Figure 5: The measured structure constant (squares) for cylinders of diameter 
N plotted against nonlinear fit for N = 4 — 8. 
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